Spin waves in disordered magnetic systems 
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Abstract. Long-wavelength spin waves in disordered magnetic systems have been investigated. 
In the framework of the Heisenberg model with magnetic dipole and exchange interactions between 
spins it is found that an additional longitudinal spin wave mode appears. This mode is characterized 
by variations of the value of the magnetic moment density. In order to analyse influence of the 
magnetic disorder on spin wave dispersion relations, the special case of volume and surface spin 
waves in the Damon-Eshbach (DE) geometry in films with magnetic disorder is considered. It 
is revealed strong influence of the magnetic disorder on surface spin waves, which consist of two 
branches - the DE mode and the longitudinal spin wave mode. Decrease of the ordering parameters 
leads to a decrease of the initial frequency and the curve slope of the DE mode dispersion curve and 
to an increase of the dispersion curve slope of the longitudinal mode. If the spin noncoUinearity is 
high, then the DE mode dispersion assumes a curve with the backward character. It is found that 
the dispersion relation of the longitudinal mode is weak temperature dependent. The developed 
model can explain the observed double-peak structure of FMR spectra in magnetic nanocomposites. 



^ ,1. Introduction 

m ; 

^ 'Amorphous magnetic materials and nanocomposites consisted of magnetic nanoparticles in 
Tj- 'insulating matrix are important types of disordered magnetic systems. Due to the amorphous 
^ structure, amorphous magnetic materials possess unusual properties such as the influence of 
QQ 'fluctuations of the exchange interaction on the magnetic state, the appearance of the roton spin 
wave branch and the giant Ai?-effect IH E], [3l H]. Nanocomposites consisting of ferromagnetic 
'metal nanoparticles embedded in an insulating matrix exhibit interesting properties such as 
^ 'giant magnetoresistance due to spin-dependent tunneling [5], El El [HI [HI [TOl [H], the occurrence 
'of additional modes in the ferromagnetic resonance (FMR) spectrum in a narrow range near 
the percolation threshold fl2[ [T3l HI], sharp increase in the FMR linewidth with decreasing 
metal content [121 [El [IS] , excellent soft magnetic properties governed by the interaction between 
ferromagnetic nanoparticles [161 [13 and the influence of inter-cluster dipolar interactions on the 
magnetic relaxation of nanoclusters [T8l [19]. The effect of the injection magnetoresistance has 
been observed on heterostructures of silicon dioxide with cobalt nanoparticles grown on GaAs 
substrates [20l[2T]. Owing to this, these heterostructures are of great interest for spintronic devices. 

In order to analyze the above-mentioned properties and effects, it is utterly important to get 
information about magnetic structures of disordered systems. One of the effective tools, which 
can be used for the characterization of magnetic structures is spin waves propagating in these 
disordered magnetics. Spin waves are very sensitive to the magnetic disorder. This exhibits in 
changes of spin wave dispersion curves. For investigation of changes of dispersion curves, one can 
use spin waves generated by antennae [22l [23l [21] or thermally excited spin waves observed by the 
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Brillouin light scattering technique [25| [261, [271 [281, [29] . In order to interpret experimental results, it 
is need to determine theoretical dependencies between spin wave dispersion curves and disordered 
parameters. This leads us to develop a model of spin waves in disordered systems. 

In this paper, in the framework of the Heisenberg model with magnetic dipole and exchange 
interactions between spins we consider long-wavelength spin waves in disordered magnetic systems. 
Averaging pseudo differential Landau-Lifshitz equations |30], in section 2 we derive the equation 
describing spin waves. In order to analyse influence of the magnetic disorder on spin wave dispersion 
relations, in section 3 we consider particular solutions for the case of spin waves in the Damon- 
Eshbach (DE) geometry. We have found that in disordered magnetic systems a longitudinal spin 
wave mode appears. This mode is analogous to plasmon modes in solid states [31], [32]. In contrast 
with disordered systems, in ferromagnetically ordered ones the longitudinal mode is degenerated 
and is not observed. Taking into account the existence of the longitudinal spin wave mode in 
disordered systems, in section 4 it is supposed that the double-peak structure of FMR spectra 
observed in magnetic nanocomposites [33l [3ll [351 ES] can be explained by excitations of usual spin 
waves and the longitudinal spin wave mode. 

2. Derivation of long- wavelength spin wave equation for disordered magnetic systems 

In order to derive equation describing long-wavelength spin waves in disordered magnetics, we 
consider the Heisenberg model with magnetic dipole and exchange interactions between spins [30] . 
We will use the approximation of effective propagators and interaction lines of the spin operator 
diagram expansion. In the framework of this approximation the local dynamics of the spin system 
is described by Landau-Lifshitz equations. The ac magnetic field originated from spins is described 
by the equation for the magnetostatic potential. The pseudodifferential Landau-Lifshitz equations 
have the generalized form and take into account longitudinal temperature dependent variations of 
the variable magnetization. Averaging over a small volume 5V and using the transformation from 
the local basis to the global one, we will obtain the average value of variations of the magnetic 
moment density. Substitution for this average value in the equation for the magnetostatic potential 
gives us the spin wave equation for disordered magnetic systems for spin waves with wavelength 

A > {6vy/^ 

In order to analyze disordered spin systems and to derive spin wave equation, we choose 
the local basis (x', y', z') such that the axis Oz' is parallel to the self-consistent field H^^^ = 
H + + H^'^'P^ (figure [ID, where H is the applied magnetic field, H^^^^ is the exchange field 
and H^'^'^P^ is the dipole magnetic field originated from neighbouring spins. The exchange field 
can be written as if^^'^^(f) = (g^J^B)'^ Ifiu{r - r'){{Su{r'))), where = {+,-,z'}, I^y is 
the exchange interaction between spins, g is the Lande factor, /is is the Bohr magneton and 
{{S^{r'))) = {{{S+{^))), ((^-(f^))), {{S.ir')))} {S± = S^^±iSy^) is the statistical average spin. It 
is supposed that the summation in this relation and in all following relations is performed over all 
repeating indices. The dipole magnetic field can be written as [30] 

Hi'''\^ = g^^BV,J2T^X((s^(^))) (1) 

r' ' ' 

= f^^-/ ^VU(^.(r-')))rfV + i7(:^)(f), 

Jv ]' ~ ' \ 
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Figure 1. Transformation from the local basis {x' , y', z') to the global basis (a;, y, z). The axis OA 
lies in the surface xOy. 



where the first term is the depolarizing magnetic field of the continuum magnetic sample, H^f^if) 
is the anisotropy magnetic field, Va is the atomic volume. Without action of ac magnetic fields at 
the site r the average value of spin {{S{f))) is parallel to the field H^'^\ In this spin orientation the 
local energy minimum is occurred. We assume that in the vicinity of this local energy point the 
exchange interaction is isotropic, 2/i°| = 2/|°l = li^J = Then, if the ac magnetic field hy{f,uj) 
acts on spins, in the local basis (a;', y' , z') the Landau-Lifshitz equations give relations between 
small variations of the magnetic moment density, m^{f,uj) = gfiBS{{S,y{r,uj))) /Va, and the field 



m±(r, uj ) 



my[r,uj) 



2-iEl^M{r)h^ 



r.uj] 



M{r)K 



r,uj] 



(3) 



where h± = l/2{h.j.^ihy), oj is the frequency, 5{{Sy{r,uj))) is the variation of the statistical average 
spin {{S^{r))), 7 = gfis/h is the gyromagnetic ratio and M(r) = g^BB{p)/Va is the magnetic 
moment density. Functions 



B{p) =SBs{Sp) 
op 

are expressed by the Brillouin function Bs for the spin S: Bs{x) = (1 + 1/2S') coth[(l + 
l/2S)x] — (1/2S') coth(x/2S'). p = l3gnBHz^\'r), 13 = l/kT, k is the Boltzmann constant, T is 
the temperature. At low temperature the derivative of the Brillouin function B^^^(p) tends to and 
the longitudinal variable magnetization ruz is negligible. The Landau-Lifshitz operators have the 
form 
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+ '^^^"^^(^) f k'^exp[ik{f-r')]m±{r',uj)d^k (4) 



E^m^{f, u) = Ljm^{f, uj), (5) 

where H^""'^) = \H + H^'^^'p)\ is the value of the sum of the external maErnetic field and the 
dipole magnetic field acting on the spin S'(r), H = (27r)^/V^, a = wVa/4:ii{gfiBy is the exchange 
interaction constant and w is the coefficient in the Fourier transform of the exchange interaction 
in the local energy minimum with respect to the space variables: i^'^^{k) = exp(— ifcr) = 

I^^\0) — wk"^ . The Landau-Lifshitz equation ([3]) is considered in the linear approximation with 
respect to B^^\p) / B{p)^ therefore, we have left the main factor oj in the operator Ez © and have 
dropped out factors proportional to B^^\p) in this operator. In order to simplify notations, we 
write the Landau-Lifshitz equations ([2]), ([3]) in the local basis in the short form 

m^(f,^) = x%^^h^{f,u}), 

where Xu'^^'' is the pseudodifferential tensor operator. 

In the considered Heisenberg model the interaction between spins contains the magnetic dipole 
interaction. From the form of the magnetic dipole interaction [30] it follows that the ac magnetic 
field hy originated from neighbouring spins is magnetostatic, i.e. it is expressed in terms of 
the magnetostatic potential ip: hy = —Vuf- In the approximation of effective propagators and 
interaction lines in the global basis (x, y, z) the ac magnetic field of spin waves is described by the 
equation for the magnetostatic potential 

— A(y9(r, oj) + 47r Vimj(r, u) = 0, (6) 

where the index i is (x, y, z). 

Spin wave equation for disordered magnetic systems is derived by averaging the Landau-Lifshitz 
equations ([2]), ([3]) over random variables and by the substitution of the average value of magnetic 
moment density variations for rui in equation ([6]). In order to substitute in equation ([6]), we must 
perform the transformation from the local basis (x', y' , z') to the global basis (x, y, z) (figure [1]) and 
average the pseudodifferential tensor operator xli/?^'' over angles ^, ■?/', 0, the exchange interaction 
I^^^ and the value of the dipole mag netic field H^'^^p'^ = \H^'^'p^\. Avera ging is performed in a small 
volume 6V of the disordered system. The volume 6V ^ Va and must contain a great number of 
spins, but its size is considerably smaller then the size of the sample. The transformation U between 
the bases {x',y',z') and {x,y,z) are expressed through the Euler angles 9, ip and [37j. Without 
loss of generality, the Euler angle can be assumed equal to 0. Then, in the global basis {x,y,z) 
the average tensor operator xlj^\ which gives relations between and hj, mi{r, u) = Xij^^hj{f, uj), 
has the form 

OO OO 7r 2iT 

xSr^ = III I Ur;x%''^U,,u{I^'\H^''^^)f{d)p{i;)s\n^ (7) 

-OO 

where indices z, j are (x, y,^), indices z/, ^ are {x',y',z ') and f{9), p{iP) are 

distributions in the volume 6V. Distributions must be normalized by relations 



/OO poo 
-OO Jo 
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f{9) sine d9 



2-K 



p{ip) dip = 1. 

We suppose that the angle ip is the random variable with the distribution p{ip) = (2 



this case, taking into account equations ([2] 
transformation U versus the Euler angles [37J 



7r)-i. In 
and well-known expressions of the 



^ cosip sinip ^ 

U = — cos 6 sin ip cos 6 cos ip sin 6 
sin 6 sin ip — sin 6 cos cos 6* y 

and integrating over ■?/' in relation ([7j), we find the average tensor operator 



X 



(av) 



where 




-oo 



-GC Fi + Dr, 
y 2Fri + D{1 - 2ri) j 



F = ]p{E-^ + E-_')M{r) 



u 



G = -^{Ep' - Ez')M{r) 

+ cos^e) sin^d^ 



(= / f{9) cos 9 sine de. 
Jo 

Substituting nii = —Xij'"'^'V jip in equation we obtain the equation describing spin waves in 
disordered magnetic systems 

A^{f, uj) + 47rVaS7V,y^(f, u) = 0. (9) 

In connection with equation ([9]) it is need to make the following remarks. 

(1) Equation ([9]) describes long-wavelength spin waves with wavelength A ^ (SVY^^. Effects 
related to short-wavelength spin waves such as the appearance of the roton spin wave branch and 
the spin wave localization are not described by equation iQ and are out of our consideration. 

(2) Relaxation of spin waves can be taken into account in the one-loop approximation of the 
considered Heisenberg model, which is the next approximation after the approximation of effective 
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propagators and interaction lines [30]. Self-energy diagrams containing n loops give correction 
terms to Green functions and interaction lines proportional to where Rmt is the radius 

of the interaction between spins. For Va/ Rf^t 1 one-loop diagrams give greatest correction terms 
in comparison with self-energy diagrams with n > 1 loops. 

It is need to note that in the approximation of effective propagators and interaction lines the 
average value of the magnetic moment density is given as 

/ M(r)Cu(/(°\if^'^'^^)rf/(°)d//('^'f\ (10) 

oo Jo 

where in the magnetic moment density M^""""^ (r) the vector r is the center of the small volume 
5V. At the full randomness of the spin disorder the ordering parameter C = and, therefore, 
M(«'')(r) = 0. But, at the same time, the tensor operator x^^'"^ dH]) takes on the diagonal form and 
does not tend to zero. 



3. Spin waves in films with magnetic disorder 

In order to clarify peculiarities of long-wavelength spin waves in disordered magnetic systems, we 
consider spin waves in the DE geometry in films with magnetic disorder. In the DE geometry the 
orientation of the applied magnetic field H is parallel to the film surface, spin waves propagate 
along the axis Ox and the wavevector q is orthogonal to the field H (figure [2]) . We consider the 
case, when the spin wavelength A ^ (SVY^^ and the magnetic disorder is due to a noncollinearity 
of spins. Therefore, we suppose that in relations (JTj) and ([H]) the distribution u(I^'^\ is equal 

to 5(/*^°^ — Iq^^)6{H''^^p^ —H^^^"^), where l'^^ is the average value of the exchange interaction between 
spins and Hq'^^^^ is the average value of the dipole magnetic field acted from neighbouring spins. 
Since the depolarizing magnetic field of the continuum in-plane magnetized film is equal to zero [23] , 
in accordance with relation ([1]) the dipole magnetic field ^q'^*^-' (r) can be reduced to the anisotropy 
magnetic field HQ"'\r). Taking into account that the film is homogeneous through the thickness d, 
we will find the solution of equation ([9]) in the form 

{Aiexp{\q\y), y<0 
A2 exp{Qy) + As exp(-Q?/), <y <d (11) 

A4exp{-\q\y), y>d 

where Q is the transverse wavevector, q = 27r/A. The magnetostatic potential (p{f,uj) and the 
normal component of the ac magnetic induction must be continuous at boundaries. This gives the 
boundary conditions 



yy 



r.uj] 



[12) 



= iTixiTVMr, ^) + (1 + 47rxi70 V,^(f, .;) , 

where d is the notation of the boundaries at y = and y = d. The boundary conditions determine 
spin wave dispersion relations. In the DE geometry two types of solutions of equation iQ exist - 
volume and surface spin waves. 




Figure 2. Scheme of the film with a magnetic disorder and distributions of the magnetostatic 
potential ip{f,Lu) of volume spin waves and surface spin waves propagating in opposite directions. 



3.1. Volume spin waves 

Volume spin waves (figure [2]) are characterized by imaginary values of the transverse wavevector 
Q. Sewing ip{r,Lo) in accordance with the boundary conditions (|T2l) . we determine the coefficients 
Ai,...,A4 in the form 

= ^4 = A2 = -^3 

and get the spin wave dispersion relation 

LOn = mH + aQMiq' + IQnl'mn + a^Aiiq' + \Qn\') + ^^m]V^' + Soo^, (13) 

where the transverse wavevector Qn = mn/d^ n = 1, 2, 3, ... is the mode number, VLh = 'jH^^"'^^ = 
+ ^o""*!; = 4:7T'yM. The term 6uJn is due to the longitudinal variations of the magnetic 
moment density in equation ([3]) caused by nonzero temperature. In the linear approximation 
with respect to B^^\p) / B{p), this term is written as 

For the case of ferromagnetic films at zero temperature, the term 5uJn = 0, the parameter of ordering 
^ = 1 and relation f|T3l) coincides with the well-known relation of spin wave spectra for in-plane 
magnetized ferromagnetic films [23] . 



3.2. Surface spin waves 

Surface spin waves (figure [2]) are characterized by real values of the wavevector Q. Taking into 
account relation (IHl) and the form of the solution (fTTI) in the interior region of the film, from equation 
([9]) we find that Q = q. Sewing ip{f,uj) in accordance with the boundary conditions (fT2l) . we can 
express the coefficients Ai,A^,A4^ in the form ffTTj) versus the coefficient A2 

_ {nl - cu^)(i + a-K) + nMi^^H + Cuj) 



Ai = A2 + A3 



(14) 
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A4 = A2exp[{\q\ + q)d\ + A-sexpKlql - q)d\ 



where the parameter k = q/\q\ = ±1 defines the direction of the spin wave propagation and 



where u = 1 — exp{—2\q\d) . The term 6uj is due to the longitudinal variations of the magnetic 
moment density in equation ([3]) and is given by 



From relation f |T4|) one can see that the magnetostatic potential if{f,Lj) of surface waves 
depends on the direction of the wave propagation - for the wavevector q > the wave propagates 
along one surface and for g < the wave propagates along the other surface (figure [2]) . The spin 
wave dispersion relation f|T5|) determines two branches of surface spin waves. In order to elucidate 
infiuence of the spin noncoUinearity on these branches, we have carried out calculations of dispersion 
curves (fT5l) without the term Su. Dispersion curves are presented in figure [3] for the magnetic film 
with the saturation magnetization AtiM = 5 kOe for different parameters of ordering, ^ and (. 
The mag netic field //("^"f) = \H + ^^")|, which is the sum of the applied magnetic field and the 
anisotropy magnetic field, is equal to 2 kOe for all curves and the gyromagnetic ratio 7 is equal to 
271 ■ 2.83 MHz/Oe. Spin waves propagating in opposite directions have identical dispersion curves. 
For = C = I spins have ferromagnetic ordering. In this case, the upper branch is the DE mode [23] , 
the lower branch at the frequency F = 00/2-71 = VLh /27i = 5.66 GHz is degenerated. Usually, the 
lower branch is not taken into account and is neglected. Decrease of the ordering parameters ^ and 
( leads to a decrease of the initial frequency of the DE mode dispersion curve and to a decrease of 
its slope. At the same time, the degeneration of the lower branch is removed. For films with the 
spin noncoUinearity of high values, the DE mode dispersion assumes a curve with the backward 
character. The limiting case of the spin disorder - full randomness - corresponds to the ordering 
parameters ^ = 2/3 and C = 0- At these ordering parameters the presented curves in figure [3] can 
be regarded as the limit of surface spin wave dispersion relations. 

4. Discussion 

Consideration of dispersion curves of long-wavelength spin waves presented above is carried out 
for films in the DE geometry and exhibits two branches of surface spin waves. In consequence of 
this consideration, it is arisen the following question: what is the nature of these two spin wave 
branches, which are observed in disordered magnetic systems? It is need to note that only one of 
these branches is observed in ferromagnets. In order to explain this problem, we consider rotation 






(15) 



Suj 



riB^^^{p)nM{co'^ - njj){2Lu'^ - 2nl - ^nH^Mu) 
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0.0 0.4 0.8 1.2 1.6 2.0 

qd 

Figure 3. Influence of the spin noncollinearity on two branches of long-wavelength surface spin 
wave dispersion curves in the film with a magnetic disorder with the saturation magnetization 
AttM ^ 5 kOe at the magnetic field = 2 kOe. The upper branch is the DE mode and the 

lower branch is the longitudinal spin wave mode. The wavevector q is normalized by the thickness 
d. The parameters of ordering are: (1) - ^ = 1, C = 1; (2) - ^ = 0.9, C = 0.7; (3) - ^ = 0.75, C = 0.3; 
(4) - e = 2/3, C = 0. 



of two noncoUinear spins in the field H^^\r) in the volume 5V (figure H^). In the common case, 
due to different orientations and values of the field H^'^\r) at different points 1 and 2 of the space 
variable r, variations 6 Si, 6S2 and their phases are different 

5i(t) = Sf^^ + 6Siexp{tuJit) 

S2{t) = S^f + SS2 exp{tuj2t), (16) 

where Sf'\ S^^ are nonperturbed spins in the points 1 and 2, respectively. Averaging in relations 
(I7j) and (IHl) in the volume 6V corresponds to the summation of rotating spins. According to relations 
(fT6l) . this leads to changes in the value of the sum of spins during their rotation (figure Hb) 

|5i(t) + S2{t)\ = [\S{'^ + 5f p + (5f) ■ 5S2) eM^^2t) + (5^f ■ 5Si) expitcuit)]'/^ (17) 

and, consequently, the value of the total spin of the volume 6V can be changed. For simplicity, in 
relation f|T7|) we restrict ourself by the linear approximation with respect to |55'i(2)|/|5^i(2)|. From 
relation ( |T7I) one can conclude that in disordered systems the average magnetic moment density 
given by relation (fTOj) can be variable. This leads to the existence of spin excitations of two 
types: (1) rotations of the total spin of the volume 5V and (2) variations of the value of this total 
spin or variations of the value of the magnetic moment density M^"''"\f). The first excitations are 
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a. 




1 2 




Figure 4. (a) Rotation of spins 5*1 and 5*2 at different points 1 and 2 of the space variable in the 
volume 6V. (b) Changes in the value of the sum of spins. 

analogous to excitations in ferromagnets. The second ones can be called as longitudinal excitations. 
They are not observed in ferromagnets. Longitudinal excitations are more expressed in the case of 
full spin randomness. In this case, ordering parameters reach the values: ^ = 2/3,( = 0, ri = l/3, 
and the average tensor operator x^""""^ ([8]) has the diagonal form. Due to the pole singularity at 
uj = riff in the operator x^"-'"') , for the longitudinal spin wave (the lower branch in figure [3]) the 
magnetic moment density variations mj = —xlj^^Vjip have high values. 

Longitudinal spin excitations are analogous to plasmon excitations in solid states. Plasmons 
are collective oscillations of the electron density [311 [32] . In plasma the Coloumb interaction 
between electrons is dominant and is determined by the electrical potential. The longitudinal 
spin excitations in disordered magnetic systems are characterized by variations of the value of the 
magnetic moment density and are described by the magnetostatic potential given by equation ([9]), 
which are analogous to the electron density and to the electrical potential in plasma, respectively. 
In this sense, the longitudinal spin excitations can be called as the plasmon-like spin wave modes. 

It is need to note differences between DE and longitudinal modes. (1) According to spin wave 
dispersion relation (fT5l) for in-plane magnetized films, the frequency u of the DE mode is higher 
than the frequency of the longitudinal spin wave mode. This increase of the DE mode frequency is 
due to the nonzero value of the saturation magnetization AttM. This leads to the nonzero value of 
Qm in relation f|T5l) . In contrast with the DE mode, the frequency of the longitudinal mode is close 
to Hh determined by the magnetic field H^"^"'^\ which is the sum of the applied magnetic field and 
the anisotropy magnetic field. (2) Owing to the inequality u > Qh for the DE mode, the term 
6ijj in the dispersion relation (JTSll depends on temperature through the Brillouin function factor 
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B^^\p) / B{p) and the frequency of the DE mode decreases with temperature increasing. In contrast 
with this, since uj — > VLh for the longitudinal mode, the dispersion relation of the longitudinal mode 
slightly depends on temperature. 

FMR spectrum is the limit case of spin waves - the spin wave wavevector q is equal to 
zero. In this limit case, due to the pole singularity aX uj = in the operator x^"''"\ for the 
the longitudinal spin wave mode in disordered magnetic systems the magnetic moment density 
variations = —Xij'"^'^jV^ ^ire not zero at 'Vjip — > 0. This leads to an additional peak in 
the FMR spectrum. In magnetic nanocomposites the double-peak structure of FMR spectra is 
observed [331 EH ESI ES] . A weak secondary peak appears in FMR spectra in nanogranular films 
composed of ferromagnetic Fe nanoparticles embedded in Si02 glass matrices p3]. FMR profiles 
for nanocomposites Ni/ZnO consisted of nanosized Ni and ZnO particles shows a secondary peak 
which appears at the high field tail part of the absorption curve [35]. Two resonance peaks are 
observed in the /i" spectra for Co/ZnO nanocomposites [M]. Microwave spectra of Ni/ZnO and 
Ni/7-Fe203 nanocomposites possess a double-peak behavior in the losses over 2-16 GHz frequency 
range [36]. Taking into account that the observed effect can be caused by oxide structures formed 
in the nanoparticle interfacial regions, another explanation of the double-peak structure of FMR 
spectra can be proposed on the base of the developed model of spin waves in disordered magnetic 
systems. According to this model, one peak is the usual rotation of the total spin. In the FMR 
spectra this peak corresponds to usual spin waves with the zero wavevector, g = 0. The second peak 
is due to variations of the value of the magnetic moment density M*^"^)(r). This peak corresponds 
to the longitudinal spin wave mode at g = 0. 

5. Conclusion 

We have investigated long-wavelength spin waves in disordered magnetic systems. The results of 
the investigations performed can be summarized as follows. 

(1) In disordered magnetic systems a longitudinal spin wave mode is appeared. In ferromagnets 
this mode is degenerated and is not observed. The longitudinal spin excitation in disordered systems 
is characterized by variations of the value of the magnetic moment density. In this sense, this spin 
excitation can be called as the plasmon-like spin wave mode. 

(2) The special case of long- wavelength spin waves in the DE geometry in films with magnetic 
disorder are considered. It is found that the magnetic disorder influences volume and surface spin 
waves. Surface spin waves have two branches - the DE mode and the longitudinal spin wave mode. 
Decrease of the ordering parameters leads to a decrease of the initial frequency of the DE mode, 
to a decrease of the slope of the DE mode dispersion curve and to an increase of the dispersion 
curve slope of the longitudinal mode. For films with the spin noncoUinearity of high values, the 
DE mode dispersion assumes a curve with the backward character. It is found that the dispersion 
relation of the longitudinal mode is weak temperature dependent. 

(3) The observed double-peak structure of FMR spectra in magnetic nanocomposites can be 
explained in the framework of the developed model. 
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